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Abstract 

We have studied positive-parity states of 28 Si using antisymmetrized molecular dynamics (AMD) 
and multi-configuration mixing (MCM) with constrained variation. Applying constraints to the 
cluster distance and the quadrupole deformation of the variational calculation, we have obtained 
basis wave functions that have various structures such as a- 24 Mg and 12 C- 16 cluster structures 
as well as deformed structures. Superposing those basis wave functions, we have obtained a oblate 
ground state band, a (3 vibration band, a normal-deformed prolate band, and a superdeformed 
band. It is found that the normal-deformed and superdeformed bands contain large amounts of 
the 12 C- 16 and a- 24 Mg cluster components, respectively. The results also suggest the presence 
of two excited bands with the developed a- 24 Mg cluster structure, where the inter-cluster motion 
and the 24 Mg-cluster deformation play important roles. 

PACS numbers: 21.60.-n, 23.20.-g 



1 



I. INTRODUCTION 



Clustering plays critical roles in excited states of p-shell and very light sc?-shell nuclei 
such as 16 O and 20 Ne . In spite of the importance of clustering in the A < 20 region, its 
role in the A > 30 region has not been studied enough. In such a heavier mass region, other 
effects such as deformations are considered to become more important than in a lighter 
mass region. Recently, based on AMD calculations, there has been discussion that both 
mean-field and cluster aspects play important roles in the excited states of 32 S, 40 Ca and 

44 Ti0,Q,B,B. 

In order to understand the mean-field and cluster aspects in the heavier system, 28 Si is 
an important case, because 28 Si has a rich variety of structures in its excited states from 
view points of both clustering and deformations. In the ground state band and the low- 
lying excited band, the coexistence of the prolate and the oblate deformation has been 
studied for quite some time. The ground state band is oblately deformed, while the excited 
band built on the 0^~ state at 6.69 MeV is considered to have prolate deformation (prolate 



ND) 



7. 



Among these states, the (3 vibration of the oblate deformed ground state 



band generates the band built on the 0^~ state at 4.98 MeV^, [lo| . 
et al. have proposed a largely de 



Furthermore, Kubono 



a 12 C( °Ne, a) reaction 
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brmed band called the "excited prolate" band based on 
13] . This band assignment, however, is not confirmed yet, 
because intra-band electromagnetic transitions have not been observed. 

In the highly excited states of 28 Si, the cluster aspects have been discussed. The excited 
states from E x = 18 to 30 MeV, which have been observed by 24 Mg( 6 Li, d), 24 Mg(a, a) and 
24 Mg(a, 7) reactions are suggested as candidates for the rotational band members that have 



an a- Mg cluster structure 
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151 ]. though a detailed theoretical study of a- 24 Mg cluster 
states has not been done. Another cluster feature of 28 Si is 12 C- 16 clustering, which has 
been intensively investigated experimentally and theoretically. Around the excitation energy 



region of 30-50 MeV, the 12 C + 16 molecular resonances have been experimentally o 



Dserved 



by elastic, inelastic, other exit channels and fusion cross sections|16l. Il7l . Il8l . Il9l . l20l . l21l . [22]. 
On the theoretical side, the 12 C- 16 molecular resonances and their relation to the low- 
lying prolate deformed states have been studied by microscopic and macroscopic cluster 



24 



25j. By the 12 C- 16 potential model, the prolate ND band and observed 



models |23l 

12 C- 16 molecular resonances are reproduced 25]. 
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We aim in the present work to investigate the nature of excited states of 28 Si, focusing 
on clustering and deformations in a unified manner. Shape coexistence and /3-vibration are 
studied, and those states are discussed in relation to cluster components of a- 24 Mg and 12 C- 
16 clustering. We also discuss the possible existence of the super deformation and developed 
o;- 24 Mg cluster states. 

In this study, we apply a theoretical framework of deformed-basis antisymmetrized 
molecular dynamics (deformed-basis AMD) + multi-configuration mixing (MCM) with con- 
strained energy variation. The deformed-basis AMD wave function enables us to describe 
both clustering and deformation phenomena in a unified manner 

as 

■ l271|. In order to study 



excited states, we first perform energy variation under two kinds of constraints to obtain 
basis functions. We shall call these constraints (3 and d constraints that are imposed on the 
quadrupole deformation (3 and distance d between the centers of mass of clusters, respec- 
tively. Then we carry out the MCM by superposing the basis wave functions to obtain energy 
levels and wave functions of the ground and excited states. The method of deformed-basis 
AMD + MCM with the (3 and d constraints has been applied already to 40 Ca and proved 
to be efficient for describing various cluster states and deformed states as well as clustering 
correlation in the deformed states. {^J In the present study of 28 Si, we naturally adopt the 
a- 24 Mg and 12 C- 16 clustering for the d constraint, because the a- 24 Mg and 12 C- 16 clus- 
tering correlations are expected to play important roles as mentioned before. The obtained 
wave functions are analyzed in order to investigate the nature of the states. 

The present paper is organized as follows: In Sec. HIl the framework of deformed-basis 
AMD + MCM is explained briefly. Results of energy variation imposing two kinds of con- 
straints and analysis of the obtained wave functions are presented in Sec. IHII In Sec. IIVt 
results of MCM and structures of low-lying states are discussed. Clustering aspects are 
discussed in Sec. |V] Finally, we give a summary and conclusions in Sec. I VI I 



II. FRAMEWORK 



We have used the theoretical framework of de 
constraints. 26] The details are presented in Refs. [3, 27 



ormed-basis AMD 
and 28. 



MCM with 
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A. Wave Function and Hamiltonian 



The deformed-basis AMD wave function is a Slater determinant of triaxially deformed 
Gaussian wave packets, 

|$int) = AVii y?2,-- - ,<Pa), (la) 

= Xu Ti), (lb) 



w*> = n r-v 



a=x,y,z 



exp 



V„\T„— 



y \ 2 



(lc) 

IXi) = ^1 T)+AI I), (Id) 
|t;) = |p) or |n). (le) 

Here, the complex parameters Zj, which represent the centroids of the Gaussian wave packets 
in phase space, take independent values for each single-particle wave function. The width 
parameters u x , i/ y , and v z are real parameters and take independent values for each of the 
x-, y- and ^-directions, but are common for all nucleons. The spin part \xi) is parameterized 
by CKj and (3i and the isospin part |rj) is fixed as \p) (proton) or |n) (neutron). The values 
{Zj, aii, f3i} (i = 1, ■ • • , A), u x , v y and v z are variational parameters and are optimized by 
energy variation as explained below. 

The trial wave function in the energy variation with constraints is a parity-projected wave 
function, 

I*') = 1± ^int>, (2) 

where ir is parity and P r is the parity operator. In this study, we will discuss positive-parity 
states. 

The Hamiltonian is, 

H = K + % + V c -K G , (3) 

where K and Kq are the kinetic energy and the energy of the center of mass motion, respec- 
tively, and Vn is the effective nucleon-nucleon interaction. We have used the Gogny D1S 



force[29|, which is one of the widely used effective forces for (beyond-)mean-field approaches. 
It is consists of the finite-range and zero-range density dependent two-body central terms 
and the zero-range two-body spin-orbit term. The form of the Gogny D1S force is given as, 
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V^n = £*g. (4) 

i<j 
2 

^ e -(fi- 2 ) 2 /ri(^ + fl^p* _ Hn pr _ M n P a P T ) 



u 12 — 

n=l 

+iWo(cri + cr 2 )k x 6{h + r 2 )k + * 3 (1 + P^h - h)p 1/3 (^p) ■ ( 5 ) 

Where the P CT and P T are exchange operators of spin and isospin parts, respectively, the <x 
is the spin operator, and the k is the operator of the relative momentum k = (px — p 2 )/2h. 
Force parameters of the Gogny D1S force are fii = 0.7 fm, W\ = —1720.30 MeV, B\ = 
1300.00 MeV, H x = -1813.53 MeV, M x = 1397.60 MeV, fi 2 = 1.2 fm, W 2 = 103.64 MeV, 
B 2 = -163.48 MeV, H 2 = 162.81 MeV, M 2 = -223.93 MeV, W = -130 MeV fm 5 and 
£3 = 1390.60 MeV fm 4 . This force has been proved to reproduce the binding energy in the 
wide mass region systematically. The Coulomb force Vq is approximated by a sum of seven 
Gaussian functions. 



B. Energy Variation 

We have performed energy variation and optimized the variational parameters included 
in the trial wave function [Eqs. ([T])] to find the state that minimizes the energy of the system 

= ^ ' , J* ; + V cnst . (6) 

Here, we add the constraint potential V cnst to the expectation value of Hamiltonian H in 
order to obtain energy-minimum states under the optional constraint condition. In this 
study, we employ two types of constraints, one on the quadrupole deformation parameter j3 
{(3 constraint) and other on the distance between clusters' centers of mass d (d constraint) 
by using the potential Vc nst , 

, vLtiP ~ P) 2 for (3 constraint, 

vcnst — \ ~ ■ {I J 

^cnst (dc m -c n ~ d Cm -c n ) 2 for d constraint. 
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Here (3 is the matter quadrupole deformation parameter, and c?c m -c„ is the distance between 
the clusters' centers of mass C m and C n , 

dc m -c n = |Rc m -c„| , (8) 
R-c m -c n = R-c m - Rc„, (9) 
Rc n , = ^-E,ec n ^ (10) 

where Aq u is the mass number of cluster C n and the expression i G C n means that the ith 
nucleon is contained in cluster C n . It should be noted that the a (= x,y,z) component of 
the spatial center of the single-particle wave function \(pi) is R ^ fT . When sufficiently large 
values are chosen for t>^ lst and v^ nst , the resultant values (3 and c?c m -c„ of energy variation 
become (3 and (i Cm _ Cn , respectively We constrain the d a -24 Mg and <ii2 C _i6 values for the d 
constraint. In each calculation of energy variation, we constrain one of the values (3, d a _24 Mg 
and <ii2 C _i6 . Other quantities such as triaxiality 7 are not constrained and optimized by 
the energy variation. 

The energy variation with the AMD wave function is carried out using the frictional 



301 ] The time evolution equation for the complex parameters Zj,aj and /3j 



cooling method, 
is 

-^ = -/.x^,(. = l,2,-,A), (ID 

where Xj is Zj, ccj or /3i, and the time evolution equation for the real parameters u x , u y , and 
v z is 

dv a dE* 

-dt=-^' {° = W)- (12) 
The quantities /ix and //^ are arbitrary positive real numbers. The energy of the system de- 
creases as time progresses, and after a sufficient number of time steps, we obtain a minimum 
energy state under the condition satisfying the given constraint. 



C. Angular Momentum Projection and Multi-configuration Mixing 

After performing the constraint energy variation for \^), we superpose the optimized 
wave functions by employing the quadrupole deformation parameter (3 and the distances 
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between the centers of mass among clusters c?c m -c„ for the C m -C n configurations, 

= (£/&l*f> 

X \ i 

+ f d K m '° n |$ dCm ~ Cn ) J (13) 

where P^k * s the P ar hy and total angular momentum projection operator, and |$f) and 
|$^ Cm_Cn ) are optimized wave functions with j3 and dc m -c„ constraints for the constrained 
values and d^ _ c respectively. The integrals over the three Euler angles in the total 
angular momentum projection operator P(, IK are evaluated by numerical integration. The 
mesh widths in numerical integration are 2n/9, n/257 and 2tt/9 for a, (3 and 7, respectively. 
Here the body-fixed x-, y- and z-axis are chosen as (x 2 ) < (y 2 ) < (z 2 ) for 7 < 30° wave 
functions and (x 2 ) > (y 2 ) > (z 2 ) for 7 > 30° ones in the case of f3 constrained wave functions. 
In the case of c?-constrained wave functions, the z-axis is chosen as the vector Rc m -c n 1 which 
connects the C m and C n clusters. The coefficients f]? K and /^ m_c " are determined by the 
Hill- Wheeler equation, 

5(«|ff|O-e«K>)=0. (14) 

Then we get the energy spectra and the corresponding wave functions that expressed by the 
superposition of the optimum wave functions, {|$f)}, {1^"" Ms )}, and {|$j 12c " lf> °)}. 



D. Single-Particle Orbit and Squared Overlap 

In this subsection, we give the definitions of single-particle orbits and the squared overlap. 
These values are useful in analysis of the calculated wave functions. 



1. Single- Particle orbits 

the mean-field description, we 



In order to analyze a AMD wave function |$i n t) from t 
have calculated Hartree-Fock-type single-particle orbits 28] given by superposition of 
single-particle wave functions \ipi) of the AMD wave function as follows. 

First, ortho normalized wave functions |</^) are obtained by the linear combination of 
\ipi). Next, the |<^) are transformed to \(pi) by the unitary transformation to diagonalize the 
mean-field Hamiltonian matrix 
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hij = + J>y fc |(« N + *°)(Wk) - ( 15 ) 

k 

where t is the 1-body operator of the kinetic energy, and the ?) N and the v c are the 2-body 
operators of the nuclear effective interaction vf 2 in Eq. (jSJ) and that of the Coulomb force, 
respectively Thus obtained \(f>i) satisfy the following equations, 



{(f>i\(f>j) = Sij, (16) 
(<Pi\i\<Pi) + Efc(^^fcl(^ N + v C ){\0j0k) ~ \0k0j)) = eiSij. (17) 

2. Squared overlap 

In order to analyze the contributions to an MCM wave function |$^) of a subset 
X = {|$f°>} (i = 1,2, • • •) of the total set of basis wave functions, squared overlap Sx is 
calculated as 

Sx=j2\(^m x) }\ 2 , (is) 

i 

where \<§>\ ) is an orthonormalized set obtained by the linear combination of \$>\ ). 

III. STRUCTURES OBTAINED BY CONSTRAINED ENERGY VARIATION 

By energy variation under the constraints on quadrupole deformation parameter (3 and 
inter-cluster distance <i Q -24 Mg and <ii2 C _i6 for a- 24 Mg and 12 C- 16 cluster structures, respec- 
tively, energy curves as functions of /3, d a .2i Mg and c?i2 C _i6 are obtained. On the curves, 
various structures appear. 

Figure [T] shows the /3-energy curves for the positive-parity states before and after the 
angular momentum projection to J n = + states. The obtained wave functions always 
have axially symmetric shapes, though the mass quadrupole deformation parameter 7 for 
triaxiality is not constrained and optimized by the energy variation. In the small deformed 
region (3 < 0.5, the system is oblately deformed and the surface has a shallow minimum. 
In the large deformed region (3 > 0.4, the system has prolate deformation and two local 
minima around j3 = 0.5 and 0.7 that we call ND and SD minima in the following discussion. 
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FIG. 1: /3-energy curves projected to positive-parity (solid lines) and J n = + states (dashed 
lines). In smaller and larger j3 regions, oblate and prolate shapes are obtained, respectively. (See 
text) 



The Skyrme SLy7 force gives similar results 



3l| . Figures Eta), (b), and (c) show density 



distributions of the oblate, ND and SD minima obtained with (3 constraint, which have no 
remarkable cluster structures. 

In order to study the change of the intrinsic structures as a functions of j3, single-particle 
orbits are investigated. The single-particle energies e, for the neutron orbits are shown in 
Fig. [3j The orbits for protons are similar to those for neutrons qualitatively, but shift by 
approximately 5 MeV due to the Coulomb energies. Two single-particle orbits are always 
degenerate due to the time reversal symmetry, and they are approximately the eigenstates 
parity except for the transitional region from ND to SD minima around (3 = 0.58. Positive- 
and negative-parity orbits are represented by solid and dotted lines, respectively, in Fig. [31 
The parity of each single-particle orbits shows that oblate and ND states has (sd) 12 config- 
urations. In the SD region {(3 ~ 0.8), negative-parity orbits intrude, and the system has the 
Ap4h configuration in which four nucleons are promoted into the p/-shell across an N = 20 
shell gap. The promotion of nucleons into the pf -shell is confirmed by the density distri- 
bution of the highest occupied single-particle orbitals at ND and SD minima (Fig. HI). The 
density distribution at the SD minimum [Fig. H](b)] is well deformed and has three nodes, 
thereby showing its p/-shell nature, while at the ND minimum, it shows a sci-shell nature. 
Therefore, the SD minimum appears as a result of the crossing of the s<i-orbits and the 
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FIG. 2: Density distribution are shown around (a) oblate, (b) ND {(3 = 0.46) and (c) SD local 
minima (f3 = 0.79) for /3-constrained wave functions, (d) type-T and (e) type-A wave functions 
for <i Q _24 M g-constrained wave functions (d a _24 Mg = 4.5 fm), and (f) type-T and (g) type-A wave 
functions for <ii2c_i6c>-constrained wave functions (cIi2q.igq = 6.0 fm). Density distributions are 
projected onto the yz-plane, where the z- and x-axis are major and minor axes, respectively. 
Symbols "+" indicate centroids of wave packets. 

p /-orbit caused by the strong deformation of the system. 

We next discuss the results obtained by energy variation while imposing the d constraint. 
Figure [5] shows the energy curve obtained by the d a _24 Mg constraint. The energy of the 
positive-parity states is shown by solid lines as functions of the inter-cluster distance. The 
energy curves before and after angular-momentum projection to J 71 = + and 2 + are given by 
the dashed and dotted lines. Two types of shapes are obtained by energy variation. One is 
the triaxial shape (denoted as type-T) [Fig. [21(d)] and the other is the axial symmetric shape 
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FIG. 3: Energies of neutrons' single-particle orbits as functions of quadrupole deformation pa- 
rameter (3. Solid and dashed lines show positive- and negative parity states, respectively. The (5 
values for oblate shapes are defined as negative values. 
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2 fm 
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FIG. 4: Density distributions of the highest-energy single-particle orbits in (a) (3 = 0.46 and 
(b) f3 = 0.79 wave functions are shown. Solid and dashed lines show density distributions of 
single-particle orbits and total wave functions, respectively. 



(denoted as type-A) [Fig. [2£e)]. The 24 Mg cluster in the a- 24 Mg cluster structures deforms 
prolately. In case of type-A, the a cluster is located on the symmetric axis of the deformed 
24 Mg cluster, while in the case of type-T, the orientation of the deformed 24 Mg cluster is 
transverse and the longitudinal axis of the 24 Mg cluster is perpendicular to the inter-cluster 
direction. The type-T wave functions are obtained in small <i a _24 Mg region. With the increase 
of <i Q _24 Mg , the structure changes from the type-T into the type-A structure. Because of the 
triaxiality of type-T wave function, two J n = 2 + states are obtained by i^-mixing as shown 
in Fig. Efb). The energy of the second J n = 2 + state of type-T and that of the J n = 2 + state 
of type-A are almost the same at d a _24 Mg ~ 4 fm, but the overlap of those wave functions is 
quite small, and hence the type-T and -A wave functions are not mixed in the result from 
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FIG. 5: d a _24 Mg -energy curves projected to (a) positive-parity and J w = + and (b) positive-parity 
and J 71 " = 2 + states. Filled circles are obtained by hand. (See text) 

MCM. As mentioned above, in the large d a .24 Mg region, the type-A wave function is favored 
and the type-T wave functions are not obtained by the energy variation. As shown later, the 
obtained type-T wave functions are found to play an important role for an a-cluster band. In 
order to check the behavior of the type-T structure in the large d a _24. Mg region and its effect 
on the band structure, we have prepared the type-T wave functions with <i a _24 Mg = 5.0 — 8.0 
fm by shifting by hand the a cluster position in the type-T wave function at <i Q _24 Mg = 4.5 
fm (filled circles in Fig. [5]). 

Figure E] shows G?i2(>i6o-energy curves projected to positive-parity and J 71 = + state. 
Two types of shapes, type-T [Fig. [2(f)] and type-A [Fig. [2(g)], are obtained due to the 
deformation of the 12 C cluster as in the case of the <i a _24 Mg constraint. Namely, in the 
type-A wave functions, the 16 is cluster located on a symmetric axis of the oblate 12 C 
cluster, while in the type-T wave functions, the symmetric axis of oblate 12 C cluster is 
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FIG. 6: di2 C .i6 -energy curves projected to positive-parity and J 71 " = + and 2 + states. 

perpendicular to the inter-cluster direction. In contrast to the G? Q _24 Mg -energy curves, the 
type-T wave functions are obtained in a small <i Q ,_24 Mg region, while the type-A structure 
is favored in a large d a -24 Mg region. The result occurs because clusters should be excited 
owing to the Pauli principle when the clusters overlap in a small d region. In order to avoid 
overlap, a symmetric axis of a prolate (oblate) cluster tends to be oriented perpendicular 
(parallel) to the inter-cluster direction in a small d region. 

In the type-A 12 C- 16 wave functions with small <ii2 C _i6 , the 16 cluster is excited and 
forms an a- 12 C-like structure, and these wave functions are similar to the /3-constrained 
wave functions at the ND local minimum. 

IV. BAND STRUCTURES 

In this section, we discuss the results obtained by the MCM calculation. 

A. MCM calculation and energy levels 

We have performed the MCM calculation by using the obtained basis wave functions. 
The adopted basis are 22 /3-constrained with j3 = 0.07-0.48 for oblate shapes and f3 = 
0.37-0.95 for prolate shapes, 15 type-T and nine type-A <i a _24 Mg -constrained wave functions 
with d a -24 Mg = 1.0-8.0 fm and 4.0-8.0 fm, respectively, and six type-T and six type-A 
<ii2 C _i6 -constrained wave functions with di2 C _i6 = 5.5-8.0 fm and <ii2 C _i6 = 3.5-6.0 fm, 
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respectively. In the MCM calculation, \K\ < 2 and ($|P^|$)/($|<3>) > 0.005 states are 
adopted as the MCM basis. Other states are omitted, because they contain numerical errors 
in the numerical integration of the angular momentum projection. The convergence of the 
MCM calculation is confirmed by changing the number of the basis wave functions. 

As shown in the energy spectra in Fig. [3, many rotational bands are constructed due to 
the coexistence of various structures. An oblate ground (g) band, a (3 vibration (vib) band, 
a ND band and a SD band are obtained. Moreover, two a cluster bands (a 2 + and a +) 
with a- 24 Mg cluster structure are found in high excited states. As for the experimental «o+ 
band, averaged values of the excitation energy are listed because those states are fragmented 
(dotted lines), 



B. Shape coexistence and j3 vibration 

The excitation energies of the ground state band and j3 vibration band members have 
good agreement with experimental data. As for the ND band, the calculated excitation 
energies are slightly higher than the experimental ones, but the calculations reproduce well 
the level spacing in the ND band as shown in Fig. [HJ which shows moments of inertia 
as functions of angular momentum. That is, both theoretical and experimental values of 
moments of inertia of the ND band are almost constant and approximately 4 h 2 /MeV. 

In order to analyze the wave functions of the ground, j3 vibration, ND, and SD bands, 
squared overlaps Sp=p t = |($ MCM |$^)| 2 with /3-constrained wave functions are used, as 
shown in Fig. M, where the $ is the value of the quadrupole deformation parameter (3 for 
the |$f ). The Sp=^ values for the band-head + states are plotted as functions of quadrupole 
deformation parameter (3. The wave functions of both the ground and (3 vibration bands 
have large amplitudes in the oblate region, and it is found that the J 71 " = 0^~ ib state ap- 
pears owing to its orthogonality to the ground state J n = 0+, which shows a (3 vibration 
mode. The ground state band amplitudes a peak at f3 = —0.4, which shows a rather large 
deformation. The ND and SD bands have large amplitudes at prolate regions {(3 ~ 0.4 and 
0.8, respectively). These results suggest the shape coexistence of the oblate and prolate 
normal-deformations and the prolate superdeformation. 
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FIG. 7: Level scheme of positive-parity states in 28 Si is shown. Left and right panels are for 
experimental and theoretical values, respectively. The experimental data are taken from Refs. 
and 132. The "ex. pr." indicates "excited prolate" band. 
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C. Superdeformed band 

The present result predicts the SD band starting form J 71 " = + state at 13.8 MeV, though 
the SD band has not been clearly identified. There are experimental works that argue for 



the existence of the "excited pro 



the state at around 10 MeV 
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ate" band with a large moment of inertia starting from 
I2I . We compare the theoretical SD band and the 

experimental "excited prolate" band. 

As shown in Figs. [1] and [SJ SD states are obtained by wave functions around the local 
minimum at (3 ~ 0.8. Reflecting the large deformation, the moments of inertia are large, and 
take a value of approximately 6 ft 2 /MeV as shown in Fig. [HJ However, the moments of inertia 
for the "excited prolate" band deduced from the band assignment of the experimental work 
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angular momentum [h] 



FIG. 8: Moments of inertia for the ground, (3 vibration, ND and SD band. Open and closed 
symbols indicate theoretical and experimental values, respectively. The definition of moment of 
inertia is J {J) = e(j+2)~-e(j) wnere ^(J) is excitation energy of the angular momentum J 
state. 
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FIG. 9: Squared overlaps with /3-constrained wave functions as functions of quadrupole deforma- 
tion parameter (3 for ground, (3 vibration, ND, and SD bands are shown. (3 values are defined as 
negative values for oblate shapes. 

are much larger than those of the theoretical SD band. The present results unsupport the 
band assignment of "excited prolate" band. In order to conclude the correspondence of the 
theoretical SD band and the "excited prolate" band, more experimental information, such 
as intra-band transitions, are required. As shown in the next subsection, large strengths for 
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TABLE I: Electric quadrupole transition strengths B{E2) are shown. Units are in Weisskopf unit 
-B(£ , 2)w.u. = 5.05e 2 fm 4 for 28 Si. Jj and Jf are the angular momentum of initial and final states, 
respectively. Experimental data and the results of the 12 C + 16 O potential model (PM) are also 
listed. The experimental data are taken from Ref. [32. 



Ji Jf experiment present PM 
intra 2+ 0+ 13.2 ± 0.3 15.0 
4+ 2+ 13.8 ± 1.3 23.2 



6+ 4+ 9.9 ± 2.5 28.6 

33.3 

2 v + ib v + ib 5.5 ± 1.3 8.7 

4b 2 vib - 14-2 

6 v + ib <ib - 15.2 

2+ D 0+ D - 41.7 45.9 

4+ D 2+ D 29 ± 5 57.4 63.9 

6 nd 4 nd > I 6 59 -l 67 -l 

8+ D 6+ D 56.3 64.9 

2 SD OgD 132.1 

4+ D 2+ D - 188.1 

6+ D 4+ D - 205.8 

8sd 6sd - 2 l 2 -6 

inter 0+ ib 2+ 8.6 ± 1.6 6.7 

2 vib 0.029 ±0.009 0.3 

2+ b 4+ 0.8 ±0.3 3.3 



the electric quadrupole transition are suggested in the present SD band. 
D. Electric quadrupole transition strengths B(E2) 

Table d shows electric quadrupole transition strengths B(E2) of intra- and inter-band 
transitions. Experimental data and theoretical values of the 12 C + 16 potential modelj^ 
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are also listed. 

Reflecting the large deformation of the ND and SD states, the B(E2) values for their 
intra-band transitions are large. The 12 C + 16 potential model gives consistent results 
with the present calculation for the intra-band transitions in the ND band. Compared to 
the experimental data, the B(E2) values are overestimated. Especially, B(E2; 6+ — > 4+) 
is much overestimated. It might be because structural changes with an increase of angular 
momentum are not represented enough in the present framework, which uses variation before 
angular momentum projection. 



V. CLUSTER CORRELATIONS 

Table HTl shows squared overlaps of the MCM wave functions with the oblate ((3ob), prolate 
ND (Psd) and SD (/3 SB ), a- 24 Mg type-T (a T ) and A (a A ), and 12 C- 16 type-T (C T ) and A 
(Ca) components. Here, /?nd and /3sd indicate subsets consisting of the /3-constrained wave 
functions with j3 = 0.37-0.58 and (3 = 0.65-0.95 wave functions, respectively. The definition 
of squared overlap is described in Sec. Ill Dl 

The values of squared overlaps of the /3ob component in the ground state band are large 
and almost unity. However, the band also contains a large amount of the «t component, 
which that means the degrees of freedom of the a- 24 Mg clustering are embedded in the 
ground state band. In the (3 vibration band, the values of squared overlaps of the /3ob 
component are also almost unity, similar to the ground state band. The values of squared 
overlaps in the (3 vibration band are also large but smaller than those in the ground state 
band, which shows that the degrees of freedom of the a- 24 Mg clustering are not activated, 
while the f3 vibration band appear by (3 vibration as mentioned in Sec. IIVBI 

The states in the ND band contains a 12 C- 16 cluster component in the present calcula- 
tion. That presence is consistent with the works of 12 C + 16 cluster models and a 12 C + 16 
potential model, where the ND band is interpreted as the lowest band of the 12 C + 16 clus- 
ter states. 25j In the present work, the 12 C- 16 higher nodal bands are not obtained, because 
the MCM basis wave functions of the type-A 12 C- 16 configuration with large <ii2 C _i6 are 
not sufficiently incorporated to describe the excitation of the 12 C- 16 relative motion. 

The SD band members contain a large a a component. Since the SD states can be 
interpreted as the 4p4h states as mentioned before, it is natural that the a-cluster correlation 
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TABLE II: Squared overlaps of obtained wave functions and oblate (Pob), prolate ND (/?nd) and 
SD (Psd), a- 24 Mg type-T (a T ) and A (a A ), and 12 C- 16 type-T (C T ) and A (C A ) components are 
shown. See text about the ND and SD. 



band J w 


Pob /3nd 


/3SD «T «A 




C A 


gs 0+ 


0.97 


0.97 






2+ 


0.97 


0.97 






4+ 


0.97 


0.93 






vib 0+ 


0.97 


0.89 






2+ 


0.95 


0.86 






4+ 


0.94 


0.53 






ND 0+ 


0.99 






0.89 


2+ 


0.99 






0.88 


4+ 


0.99 






0.88 


SD 0+ 




0.94 0.88 


0.15 




2+ 




0.94 0.88 


0.16 




4+ 




0.94 0.88 


0.16 




a + + 


0.13 


0.72 






2+ 


0.25 


0.76 






4+ 


0.05 


0.81 






a 2+ 2+ 


0.17 


0.98 






3+ 




1.00 






4+ 


0.07 


0.99 







is enhanced in this band. Therefore, the SD states are expected to be observed in a-transfer 
reactions. It should be noted that the SD state contains only the component of the type-A 
a- 24 Mg structure that has the very longitudinal structure. This property is in contrast to 
the case of the ground state band which contains only the type-T a- 24 Mg component. 

The a + and the a 2 + bands have the characteristic feature that they have large amounts 
of the type-T a- 24 Mg component, while the other components are small. Especially the 
amount of the a- 24 Mg component in the a 2 + band is quite large and nearly equals unity. 
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distance between a and 24 Mg clusters [frn] 




distance between a and 24 Mg clusters [frn] 



FIG. 10: a- 24 Mg cluster structure component S a of the (a) J T = + states in the gs and a + 
bands and (b) J 71 " = 2 + states in the g, a + and a 2 + bands as functions of distance between a and 
24 Mg clusters. 

This value indicates that this band is formed by almost pure type-T a- 24 Mg wave functions. 
That is the reason why we call these two bands the a cluster bands. For limitations of the 
present model space such as the bound state approximation, the widths or fragmentation of 
the a + are not reproduced. 

Recall that the type-T a- 24 Mg component is contained in the ground state band as 
well as the and the a<i+ bands. Therefore, it is expected that these bands might be 
understood as cluster excited states built on the ground state band. In order to analyze the 
cluster features of these bands, we here discuss the overlap of the states with the a- 24 Mg 
cluster wave functions in more detail. Figure [10] shows the type-T a- 24 Mg cluster structure 
component S a of (a) the J 71 = + and (b) 2 + states in the ground, a + and a 2 + bands as 
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functions of the distance c? a _24 Mg . As shown in the figure, the amplitudes for the J n = 0+ 
and 2+ in the ground state band are concentrated in the small d a _24. Ms region, while those 
for the J n = 0+ Q+ and 2+ + states are compressed in the small distance region compared 
to the ground state band, and they have a peak at d a _24 Mg ~ 5 fm. This figure shows the 
typical a- 24 Mg higher-nodal nature of the band built on the ground state band owing to 
the excitation of inter-cluster motion. As mentioned before, the candidates for the members 
of the a-cluster band have been observed by 24 Mg( 6 Li, d) and 24 Mg(a, a) reactions. The 
observed K n = + band may correspond to the a + band though the excitation energies are 
slightly overestimated by the present calculations. 

As for the a 2 + band, the 2„ 2+ state shows a large amplitude of S a in the small d a .24 Mg 
region, similar to the J 71 " = 2+ in the ground state band. It shows that the a 2 + band is 
regarded as a counter part of the ground state band owing to A^-mixing because of the 
triaxial deformation of the type-T a- 24 Mg cluster structure. In particular, the a 2 + band has 
a, \K\ = 2 feature of the type-T a- 24 Mg cluster structure. Here, K is defined with respect to 
the z-axis set to the inter-cluster direction of the a and 24 Mg clusters, and the orientation 
of the prolate 24 Mg cluster is perpendicular to the z-axis as mentioned before. Therefore, 
the a 2 + band is interpreted as the a- 24 Mg cluster band with the cluster core excitation of 
24 Mg(2 + ). On the experimental side, there is no established band with K n = 2 + . In order 
to search for the a 2 + band, observation of unnatural parity states might be helpful. 

28 Si has a variety of deformed bands, and those states contain cluster components. Other 
deformed states and cluster correlations such as a hyperdeformed states and a- 20 Ne-a clus- 



tering are also attractive issues [3 3l. 
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VI. SUMMARY AND CONCLUSIONS 

Positive-parity states in 28 Si have been studied using the deformed-basis AMD + MCM 
focusing on clustering and deformation. The experimental energy levels in the low-energy 
region are reproduced well by the present calculations. The oblately deformed ground state 
band and prolately deformed excited band are reproduced, and the result shows shape 
coexistence. The /? vibration band also appears because the oblately deformed state is soft 
against quadrupole deformation. A superdeformed band is suggested in the present results. 
The SD band is described by the (sc?) 8 (p/) 4 configuration. If the suggested SD band of 28 Si 
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is observed experimentally it should be the superdeformation of the lightest srf-shell nucleus. 
Existence of largely deformed band ( "excited prolate" ) has been proposed experimentally, 
however, we cannot assign the experimental band with the theoretical SD band because the 
experimental moment of inertia is not consistent with that of the calculated SD band. More 
experimental data such as electric or magnetic transitions are requested. 

The cluster bands, «o+ and «2+ bands also have been obtained. These bands contain 
significant a- 24 Mg cluster structure components. The band is regarded as the higher- 
nodal band owing to the excitation of inter-cluster motion, while the «2+ band is interpreted 
as the K w = 2 + band due to the triaxiality of the a- 24 Mg cluster structure. 

It is found that cluster components are significantly contained in the low-lying deformed 
states as well as the cluster bands. Namely, the gs and SD bands contain a significant a- 
24 Mg cluster structure component, and the ND band contains a 12 C- 16 cluster structure 
component. Those results are analogous to situations of other sti-shell nuclei such as 32 S 
for which 16 0- 16 correlations in the SD band have been suggested j^, 0], and 40 Ca where 
a- 36 Ar correlations in the ND band and 12 C- 28 Si correlations in the SD band have been 
discussed. (3, 6] 

We have shown the importance of clustering effects as well as deformation effects in 
low- lying states of 28 Si. That finding is consistent with the recent full-microscopic studies 
that suggested that both clustering and deformations play important roles in the wide-range 
s<i-shell region. 
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